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This paper describes the determination of the first four semi-invariants of the 
distribution of the average, over an arbitrary time interval, of traffic carried by a 
telephone system with an infinite number of trunks, during a period of statistical 
equilibrium. Both finite and infinite numbers of independent call sources are con- 
sidered, and the distribution function of call holding times is left general. 

1. Introduction 

FOR mathematical studies of telephone traffic, like those of call loss or 
delay which are used in trunking engineering, the traffic is considered 
as a flow of probability in time. In the period of most importance, the busy 
hour, this flow is usually regarded as stationary ; that is to say, the proba- 
bility of a given number of busy trunks, or the probability of delay of an 
incoming call (or any other probability of the system which comes in ques- 
tion) is taken as independent of the particular moment in the busy hour at 
which the system is examined. The system is said to be in statistical equilib- 
rium. 

For such theoretical studies, the statistical quantities which determine 
these probabilities, like the rate at which calls appear, are of course taken 
as given, but in the application they must be determined by observations, 
such as those being taken in the current extensive program of traffic meas- 
urements. Here a difficulty appears. To abridge the extensive amount of 
observational material, either measurements are made of traffic averages 
over periods small compared to the busy hour (but not small enough to be 
neglected) or the measurements of continuous recorders are averaged by 
hand. It may be noticed here that for application of the results given below 
the traffic averages obtained by measurements must be those of a con- 
tinuous device which records all traffic changes and not, as in some measuring 
devices, those obtained from a number of "looks" at points within the 
averaging interval. But to use these measurements in determining the 
traffic parameters by standard sampling theory, a corresponding theoretical 
study of the averages is necessary. 

Such a study, within limits to be described presently, is given here. No 
attempt is made to describe the sampling studies possible from the results 
reached. These seem to be of many kinds, not necessary to describe, but for 
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concreteness it may be mentioned that the most important, at the moment, 
seems to be that of setting confidence limits for the average traffic. 

The most important of the limits to this study are those implied by the 
assumptions of statistical equilibrium with fixed average, and an infinite 
number of trunks. The former limits application to periods in which, roughly 
speaking, average traffic is neither rising nor falling; the latter is justified 
only by the extreme mathematical difficulties produced by assuming other- 
wise. The traffic variable is the number of busy trunks in a period of statisti- 
cal equilibrium. For pure chance call input, the call holding time character- 
istic is left arbitrary throughout the development, but main interest lies in 
the two extreme cases of constant holding time and exponential holding 
time, which are examined in detail.* For calls from a limited number of 
sources, results are obtained only for exponential holding time. 

More precisely, if N(l) is the random variable for the number of busy 
trunks at time /, the variable studied, the average number of calls in an 
interval of length T, is 

M{T) = \ r f N{1) dt (1) 

l Jo 

The question is: What are the statistical properties of M(T)? 

The results given are the first four cumulants (semi-invariants) of M{T), 
which seem to have the simplest expressions. For the convenience of the 
reader it may be noticed that the first cumulant is the mean, the second 
the second moment about the mean which is the variance, the third the 
third moment about the mean, and the fourth the fourth moment about 
the mean less three times the square of the variance. 

In all cases the mean of M{T) is the mean of N{t) and for pure chance call 
input is called b, the average number of calls in unit average holding time, h. 

The other cumulants for pure chance call input, k n , have the general 
expression 

ti(n - 1) 



with 



kn = b "^ 1} f o dxg(x){T - .v).v- 2 ; n = 2,3,4 



* F. W. Rabe [6J has reported results for these two cases for relatively long averaging 
intervals, which are verified below. I owe my interest in this problem to a report on Rabe's 
work made by Messrs. Gibson, Hay ward and Seckler in a probability colloquium at Bell 
Telephone Laboratories initiated and directed by Roger Wilkinson. 
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and/(0 the probability that a call lasts at least /, that is, the distribution 
function of holding times. The specializations of this, for constant holding 
time and exponential holding time, appear in section 4. The results for 
finite source input have a similar character. 

The procedure in obtaining these is as follows. The cumulants are de- 
termined from the ordinary moments (about the origin) and the latter 
are determined by the integration of expectations. Thus the first moment, 
the mean is determined from 

E[M(T)] = % I E[N{1)] dt = E[N(t)} (2) 

1 Jo 

where E(x) is written for the expectation or mean of x. 
Similarly the second moment is given by 

E[M\T)\ = i [ ( E[N(t)N(u)] dt du (3) 

I ~ Jq Jo 

and so on for higher moments. Correlation effects appear in (3) in 
E[N(t)N(u)] and are included in the development by formulation of transi- 
tion probabilities, that is, those probabilities determining the traffic flow 
in time. The transition probability Pjk(l) is defined as the probability of 
transition in / from j calls in progress (busy trunks) to k calls in progress, 
and fixes the inter-relatedness of call probabilities at different time epochs. 
Only for large values of / are these probabilities independent. 

Hence, the first task is to determine these simple transition probabilities, 
then those of double and triple transitions, then the expected values of 
pairs, triples and quadruples of numbers of busy trunks, and finally the 
moments. 

2. Transition Probabilities 

For exponential holding time, and infinite sources, infinite trunks, these 
probabilities have already been determined by Conny Palm [5]. Palm's 
work has been summarized both by Feller [1] and by Jensen [3], and de- 
scribes the whole process, not merely the equilibrium condition. For the 
equilibrium condition, a different procedure,* similar to that used by 
Newland [4] for another purpose, allows the assumption of a more general 
holding time characteristic. 

* Thanks are due S. O. Rice for suggesting this, as well as for many corrections and 
improvements. I also have had the advantage of a careful reading of the mss. by E. L. 
Kaplan. 
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For infinite sources, and calls arriving individually and collectively at 
random with average density a, the well-known formula for the probability 
that exactly k calls arrive in time interval / is the Poisson 

7r k (t) = e-'W/kl (4) 

Then, if /%•/(<; k) is the conditional probability of transition from i to j 
when k calls arrive in time /, 

PM = EP«M)ir*«) (5) 

fc=0 

Consider -P ,■>(/; 0), that is the (conditional) transition probabilities when 
no calls arrive. Let the probability that a call lasts at least / be/(/), so that 
the average holding time h is given by 

h = [ u[-f'(u)] du = [ /(«) du (6) 

Jo Jo 

The i calls initially in process are independent of each other. Select one of 
them and suppose the time from its arrival (its age) is h ■ Then the proba- 
bility that it will also exist / units later is the conditional proba- 
bility/^ + h)/f{t\). Since in equilibrium conditions all moments of arrival 
have equal probability, the corresponding probability for an arbitrary call is 

g(t) = f f(t + h) dh - | f(h) dt = i ^ /(«) du (7) 

Hence the transitional probability P,y(/; 0) is the binomial expression 

P i] (f,0) = (*V(1 - *)«-' (8) 



and its generating function is 

Pt(t, x\ 0) = X>«fc O)*' = [1 + (x - l)gV (9) 

In (8) and (9), for brevity, the argument of g is omitted. 

Now, suppose one call arrives in interval t. The moment of arrival is 
uniformly distributed in /; that is, if «i is the moment of arrival, 

Pr(u < tii < u + du) = du/l 

and the probability that a call arriving at an arbitrary moment will be in 
existence at time / is, say, 

Q(t) = ff(t - u) * = i /"'/(«) ^ - 7 (1 " fW) (10) 

Jo I t Jo I 
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The corresponding generating function is 

1 - Q(l) + xQ(t) = 1 + (.r - 1)Q(/) 

and, since calls arriving are independent, the generating function for k 
calls arriving is 

[1 + (* - \)Q\" 

and 

PAt, x; k) = [1 + (x - l)f]ll + (x - l)Q)" (11) 

Hence, finally by (5), 

Pi{t;x) = E-PtfW*' 

= [i + c* - \) g ? £ [i + (* - 00* e ^-- 

= [1 + (x - Dgl'exp (x- l)at<2 

- [1 + (x - l)g]*'exp (x - 1) a/z (1 - g) (12) 

The last step uses (10). 

This is the generating function for the simplest transition probabilities, 
and is quite like Palm's result; indeed, for exponential holding time g = / = 
e~" h . The probabilities themselves are obtained by expansion of the generat- 
ing function in powers of x, or by substituting g for e~' h in Palm's result. 
But they are not needed here; the generating function is most apt for deter- 
mining the averages of interest, as will appear. 

Before going on to the other transition probabilities, it is interesting to 
notice certain checks of equation (12). In statistical equilibrium the traffic 
has Poisson density (Palm I.e.) that is, in the present notation 

Pr(N(t) = k) = <fV/£! 

where b = ah. This of course is independent of time. Then, if N(Q) has this 
density, so should N(t) as determined from A7(0) and the transition proba- 
bilities implicit in (12). This is verified by 

-b,i 

£ P i (t,x)e- b b i /il = exp (x - 1)6(1 - g) E [1 + (x - l)gY '— 

= exp [(x - l)b(l - g) - b + b + (x - \)bg\ (13) 
= exp (x — l)b. 
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Also, g(0) = 1 and g(») = () so that 

P,(0,.v) = [1 + (.v- 1)]'= *' (14) 

P,-(*>,.t) = exp (.v - 1)6 (15) 

showing that in zero time no transit to another state is possible, and in 
infinite time the equilibrium probabilities are reached no matter what the 
initial state has been. 

Finally, in a Markov process (cf. Feller [2], Chap. 15) the simple transi- 
tion probabilities alone are needed since 

P u ,(l, u) = PMPjkiu) 
A test for this is the Chapman-Kolomogorov equation, namely 

Pik{t + «) = E PuOPjM 

j 

Using (12), the corresponding relation of generating functions is 
[1 + (x - l) f (/ + u)Y exp b(x - 1)[1 - g(l + »)] 

= [1 + (* - l)*(0*(«)]'«pft(* - 1)[1 - g(t)g(u)}; 
so the process is Markovian only if 

g(l + «) = g(l)g(u) 
which is true only for exponential holding time. 

For the next transition probability P,jk(t; u), consider first the condition v 
in which no call arrives in the whole interval / + u. As before 



PiAD = (j)gi(i - a) 



where for convenience g t is written for g(l). For the next transit, however, 
there is a difference, namely 

. fc / \ j-k 

Pjdu) = " " '"" 



^©(tyc-ty 



since gt+u/gt is the conditional probability that a call which has lasted / 
will last u more; Pjk(u) is the conditional probability of a transit from j 
to k in u, given the transit i toj in /. 

The generating function for the double transition probabilities in this 
case is, then, 

E E P<»ft«;0)*y = [1 + (-v - l)g t + x(y - 1 W" (16) 

i k 

Now suppose a single call arrives at random in interval /. As before, the 
probability that it will occupy a trunk at time / is Q{t) = A/ -I (l — g(t)) 
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and the conditional probability that it will also occupy a trunk at time 
/ + u is 

\ ff(l + u- x)dx + Q(t) 
t Jo 

or 

M -_* + *) _ R(/iU)t ^ 

The corresponding generating function, with x and y the indicators of calls 
at / and / + u, resp. is 

1 - Q(t) + Q(W ~ W, «*)]* + Q(')R(t, u)xy 



or 



1 + (x - 1)(1 - g(t))h/t + x(y - l)\g(u) - g(t + u)]h/t 

The generating function for c calls in this interval is this expression raised 
to the c'th power, since calls arrive independently; and since c calls arrive 
with probability e~ a '(at) c /c\, the generating function for calls arriving in 
this interval is 

£ [1 + (.v - 1)Q + x(y - \)QR\ c e- at (al) c /c\ 

= exp b[(x - 1)(1 - g(t)) + x(y - l)(g(«) - g(l + «))] 

For brevity Q and R have been written for Q(t) and R(l, u). 
Finally the generating function for calls arriving in /, / + u, is 

expb(y- 1)(1 - g(u)) (18) 

Hence 

Z E Ptfk«, «)*'/ = [1 + Cv - l)g(t) + x(y - l)g(t + »)]' 

; ft 

• exp &[(* - D(l - g(0) + (y - D(l - g(u)) (19) 

+ x(y - l)(g(u) - g(l + «))] 

By similar argument, the generating function for triple transition proba- 
bilities is 

EZE Pijktit, u, v)x , y k z l 

i k i 

= [1 + (.v - \)g, + .v(y - l)g l+ll + acyCi - lW, + „]'' 

• exp b((x - 0(1 - g t ) + (y - 0(1 - g-„) + (20) 
(z - 1) (1 - gj + .v(y - l)(g„ - g t+u ) + 
y(z - l)(g, - £„+,) + xy(z - \)(g u +v - gt+u+v)) 
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3. Expected Correlations 

Correlation expectations, like E[N(t)N(u)] in equation (3), are needed 
for evaluation of the moments of M(T). They may be determined from the 
transition probability generating functions, if it is agreed, as a matter only 
of convenience, that the time epochs /, u, v, etc. are in that order (I < u < 
v< ■••)■ Since, on the assumption of statistical equilibrium, the call 
probabilities at the first epoch /, are independent of its value, as already 
noticed, this value may be taken as zero without loss of generality. 

Thus for the second moment it is sufficient to determine 

<p(u) = E[N(0)N(u)] = £ ipi T,JPa(") (21) 

with p t = Pr[N(0) = i] = «-»6'/*l 
Write 

G u (x, y) = E pt*' E PiMy 

By (12), this is the same as 

G u (x, y) = exp b[x - 1 + y - 1 + (* - l)(y - 1)*(«)1 
or 

Hu(x, y) = G u (x + 1, y + 1) = expb(x + y + xyg(u)) 
and 



<p{u) = 



a 2 F 



dxdy I>v =o (22) 

= 6" + bg(u) 
In the same way the second order correlation expectation, that is 
<p(u, v) = E[N(0)N(u)N(u + »)], 
is obtained from 

Gu,v(x, y, s) — £ frs* E Z) -P<ifc(«> »)yV 
and 

#«.»(*, y, z) = G B ,,(s + 1, y + 1, z + 1) 

= exp b(x + y + a + »yg( M ) + yzg(») + *(y + l)zg(« + »)) 
Hence 

?(«, i) = b* + Hg(tf) + g(*) + g(u + »)] + bg(u + p) (23) 



(24) 
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Finally, the third order correlation turns out to be 
<p(u, v, w) = E[N(0)N(u)N(u + v)N(u + v + w)] 

= V + b 3 [g(u) + g(v) + g(w) 

+ g(u + v) + g(v + w) + g(u + v -f- u>)] 

+ Hg(" + ») + g(* + «0 + 2g(« + d + w)] 

+ b-[g(u)g(w) + g(« + v)g(v + w) 

+ #(*>)£(« + v + w)] + 6g(« + v + w) 

As will appear, the arrangement of terms in (22), (23) and (24) corresponds 
to the expansion of ordinary moments in terms of cumulants (semi-in- 
variants); e.g. (24) corresponds to *» 4 = b 4 + 6b 2 k 2 + 4bk 3 + 3kl + £ 4 
with &,- the i'th cumulant (for the Poisson of mean b, ki = b). 

4. Moments 

Moments are obtained from these results by integrations. As already 

noted, equation (2), the first moment is b for any holding time distribution. 

Since there are two ways of ordering the epochs /, u, the second moment is 

E[M\T)] = lj dtj du<p{l - u) 

= * 2 +^o [ dt f dugQ-u) (25) 

1 - Jo Jo 

= b 2 + ^j* dxg{x){T -x) 

The last step is by the formula for reversing the order of integration indi- 
cated by 



/ dt I du = I du I dt 

Jo Jo Jo Ju 



The variance or second central moment, which is also the second cumulant 
k 2 , is then 



Var [M(T)] = E[(M(T) - bf] 
= E[M\T)] - b 2 



^^j'dxgixHT-x) 



(26) 
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Since there are 3 ! = 6 ways of ordering 3 epochs, the third moment may 
be written 

E[M\T)\ = J, f dl f du [ dv <p(t - u, u - v) 
I 6 Jo Jo Jo 

= b 3 -{-~ f dt fdu r di\ g {t - w) + g(« - v) + gn - v)\ 

T 6 Jo Jo Jo 

+ % f dl fdu f dvg(t- v) 

1 8 Jo Jo Jo 

Here the first triple integral is immediately evaluated by use of the identity 

2 f dt I du I dv[g(t - u) + g(« - v) + g(t ~ v)] 

Jo Jo Jo 

= f J [ dtdudv g(\ t — v\) 
Jo Jo Jo 

= 2T I dxg(x)(T-x) 

Jo 



T 

,3 



= r h/b 

The last triple integral, by successive inversions of integration order, turns 
out to be 

~ z [ dx g{x){T - x)x 
Hence finally 

E[M\T)] = b 3 + 34*, + % f dx g(x)(T - x)x (27) 

1 * Jo 

and 

k 3 = E[(M(T) - b) z ] 

= E[M*(T)] - Zb E[M\T)\ + 2b z 

= E[M\T)} - 3b k 2 - b z < 28 > 



= 6 j 3 j o T dxg(x)(T-x)x 
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The fourth moment is given by 

E[M\T)\ - ~. I dt I du f dv [ dwcpQ - u, u - v, v - w) 

I 4 Jo Jo Jo Jo 

= b 4 

+ P j^ 3 £ f a f" [ dt du dv dw\g(t - u) + git - v) 

+ gQ - w) + g{u - v) + g{u - w) + g(v - w)] 

»T ft /.it fV 

+ b 2 / / / dt du dv dw [g(t - v) + g{u - w) + 2g(t - w)] 

Jo Jo Jo Jo 

fT ft fU fV 

+ b 2 I II I dt du dv dw[g(t — u) g(v — w) + 
Jo Jo Jo Jo 

g(t - v)g(u — w) + g(t — w)g{u — v)] 

+ b I I I I dldudvdw[g(t - w)] 
Jo Jo Jo Jo 

Employing the identities 

• T nT r T P T 



4 / / / dtdudv dw [g(t — u) + g(t - v) + g{t — w) 

Jo Jo Jo Jo 

+ g(u - v) + g{ll — w) + g(v — w)\ 

fT ft fU rv 

= I I I I dt du dv dw g{\ I — u |) 
Jo Jo Jo Jo 

= IT 2 j T dx g(x)(T - x) = T A h/b, 

Jo 

/»7* pt p u pv 

8 / / / I dt du dv dw [g(t — u) g(v — w) + g{t — v) g(u — w) 
Jo Jo Jo Jo 

+ g(t - w) g{u - v)] 

= / / / I dt dw dv dw g(\ I — u \)g{\ v — w\) 
Jo Jo Jo Jo 

= tkl/b 2 , 



= 4 / dx g(x)(T - v) 



and successive inversion of order of integration, the final result turns out 
to be 
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E[M*(T)] = b* + Gb 2 k 2 + 46* 3 + Ski 

+ ™[dxg{x){T-x)x 2 



(29) 



and 



(30) 



h = E[(M(T) - b) 4 ] - 3E[(M(T) - b) 2 \ 

It is a tempting surmise that 

kn = b n{n ~ l) f dx g(x)(T - x) x n ~ 2 
I " Jo 

but this has not been proved. Note that for g(x) = 1, k n = b, the cumulant 
of the Poisson, as it should. 

For the two cases of chief interest, constant and exponential holding 
times, the function g(x), in average holding time units (that is, x = t/h) is 
given by 

c.h.t. g(x) = 1 - x x < 1 

= x > 1 

e.h.t. g(x) = e~ x 

and the results are as follows: 

Constant Holding Time 
Cumulant T < 1 T > 1 

k 2 b(l - T/3) bT-'i). - 1/3T) 

h b(l - T/2) bT-^1 - 1/2T) 

h b(l - 3T/5) bT~ 3 (l - 3/5T) 

Exponential Holding Time 
k 2 2bT~ 2 [T - 1 + e~T] 

h 6bT~ 3 [T - 2 + (T + 2)e-r] 

h UbT-^lT - 6 + (T 2 + 4T + 6)e~T] 

It may be worth noting that, if the surmise is correct, for constant holding 

time 



k n = b 



1 - ^-—4 T I T I 

n + 



b 



n -1 1] 
n + 1 T] 



1 -'^-rA4,\ T> 1 
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Fig. 1. — Comparison of variances of average traffic for constant and exponential hold- 
ing times. 
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Fig. 2. — Cumulants k 2 , k3, and Id for constant holding time. 
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Fig. 3. — Cumulants k 2 , k3 and k* for exponential holding time. 
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and for exponential holding time 

kn = b n(n ~ 1} [(« - 2)! T - (» - 1)! + e~ T (T + a)"" 2 ] 
where in the last term (T + a)" -2 is a symbolic expression or shorthand for 

,rr, , \n-2 V* I 11 ~ 2\ ™n-2-m 

and a m = (w + 1)!; e.g. 

(r + a y = r 3 + 6P + isr + 24 

For small values of T, the two cases coalesce (er x tt 1 — x) and at T = 
approach b as they should. For large values of T, and constant holding 
time, 

K~h/T~* t (»= 2,3,4); 

for exponential holding time 

k n ~ n\b/T"-\ (m= 2,3,4). 

For w = 2, these results agree with Rabe [6]. 

As r increases, for either holding time, the cumulants are progressively 
smaller, and the approximation of the distribution of M{T) by a normal 
curve (which has all cumulants, except the first and second, zero) improves. 
This is what follows from the central limit theorem if the subdivision of T 
into a large number of intervals results in mutually independent random 
variables (cf. Rice [7] 3.9). 

Figure 1 shows a comparison of the variances (k 2 ) for the two holding 
time cases. Figure 2 shows a comparison of the cumulants ki , kz and ki for 
constant holding time, and Fig. 3 shows the same thing for exponential 
holding time. 

5. Finite Sources — Exponential Holding Time 
The generating function for transitional probabilities for N subscribers, 
each originating calls independently with probability X, and for exponential 
holding time, as given by Jensen (I.e.) is as follows: 

P& x) = [1 + qi (x - 1)]'[1 + q (x - l)]»-» (31) 

with 

qo = p — pe 

qi = p + q " 

p = i_ g = x/(X+ 7 ) 

7 = 1/7; 
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It should be noticed that for / = «> , q = q x = p and 

P,(oo,*) - [1 + p(x- 1)]» (32) 

The right hand side is the binomial generating function and, as independent 
of i, is the generating function for the statistical equilibrium probabilities; 
that is 



Pr [N(t) = k] = 



ffi>v- 



Also the process is Markovian since 
Z x" 52 PiMPjM = Z>«0)[1 + gi„(* - I)] 7 '!! + g 0u (x - i)] N - j 

h i i 

= [1 + (qou + QuQiu — quqou)(x — 1)]' 
[1 + (qou + qotqiu — qotqou)(x — l)] N ~ 
and 

qcu + quqiu — qitqou = Qi,t+u 

?0u + qotqiu — qotqou = qo,t+u 

Here it has been convenient to indicate by the double subscript the de- 
pendence of qo and q\ on a time variable. 

Moments are obtained by the process given in detail for the infinite 
source case. For brevity it is convenient to use the binomial cumulants 
which are as follows 

K2 = Nfq 

k 3 = Npq{q - p) 

k 4 = Npq(l - 6pq) 
and the modified time variable T\= (X + y)T. Then the results are 
k 2 = 277 2 k 2 [7\ - 1 + e~ Tl ] 
h = 67T 8 /cs[ri - 2 + (Ti + 2)e~ Tl ] 
k, = 12T?((K4 + k11T 1 )[2Ti - 6 + (Tl + 47\ + 6)e~ Tl ) 

- kIN~\\- (T\+2)e-^+e-^\) 

These of course bear a strong resemblance to the infinite source case (ex- 
ponential holding time), to which they converge. 
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